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ABSTRACT 

Significant  effects  of  boundary  layer  characteristics  on 
cavitation  and  the  effect  of  unsteady  cavitation  on  noise  have 
been  widely  observed  experimentally.  In  order  to  better 
understand  these  effects,  an  experiment  with  a  two-dimen¬ 
sional  hydrofoil,  under  sinusoidal  pitching  oscillation,  was  con¬ 
ducted  in  the  36-inch  water  tunnel  at  DTNSRDC .  Three  pm 
hole-type  pressure  transducers  were  installed  on  the  foil  sur¬ 
face  to  measure  pressure  fluctuations  and  a  hydrophone  was 
used  to  measure  the  radiated  cavitation  noise.  Two  subjects 
are  discussed  in  this  paper:  l)  relationship  of  boundary  layer 
transition  and  pressure  fluctuations  with  the  hydrofoil  in 
oscillation,  and  2)  noise  generated  by  unsteady  cavitation. 

The  magnitudes  of  pressure  fluctuations  in  transition 
and  turbulent  regions  are  found  to  be  independent  of  oscilla¬ 
tion  amplitude  and  frequency.  However,  the  development  of 
boundary  layer  and  occurrence  of  transition  are  delayed  with 
an  increase  in  oscillation  frequency.  The  influence  of  hydrofoil 
oscillation  on  cavitation  inception  is  also  discussed.  With  the 
occurrence  of  leading  edge  sheet  cavitation,  a  significant  ef¬ 
fect  of  foil  oscillation  on  cavitation  noise  is  measured.  At  high 
reduced  frequencies,  extensive  cloud  cavitation  is  formed  dur¬ 
ing  the  final  phase  of  sheet  cavity  collapse  and  the  noise  level 
is  significantly  increased. 

PREFACE 

An  earlier  paper  by  the  present  authors  at  the  12th 
ONR  symposium  dealt  with  two  subjects,  the  inception  of 
cavitation  on  a  two-dimensional  foil  and  the  physics  of 
leading-edge  sheet  cavity  stability  and  subsequent  formation 
of  cloud  cavitation.  The  method  developed  to  predict  incep¬ 
tion  was  based  on  Giesing's  unsteady  airfoil  theory  and 
assumed  that  viscous  effects  were  of  secondary  importance 
The  present  paper  will  explore  the  importance  of  the  bound¬ 
ary  layer  development.  In  addition,  as  a  continuation  of  the 
early  study  of  cavitation  instability,  the  importance  of  cloud 
cavitation  in  producing  acoustic  noise  will  be  presented. 

1.  INTRODUCTION 

Cavitation  frequently  is  unavoidable  on  the  lifting  sur¬ 
faces  of  hydrofoils  and  marine  propellers.  For  many  years, 
scientists  and  engineers  have  dealt  with  the  problems  of 
cavitation-induced  noise,  structural  vibration  and  erosion 
associated  with  the  operation  of  marine  vehicles  and  hydraulic 
equipment.  All  three  problems  are  related  to  the  inception, 
growth,  and  collapse  of  vapor  cavities.1 

Inception  of  cavitation  in  a  fluid  is  the  condition  under 
which  cavitation  is  first  detected,  either  visually  or  acoustical¬ 
ly  It  had  often  been  assumed  that  when  making  analytical 
predictions,  cavitation  inception  occurs  immediately  after  the 
static  pressure  on  the  body  surface  becomes  equal  to  or  less 
than  the  vapor  pressure  of  the  fluid.  A  research  model,  the 
International  Towing  Tank  Conference  (ITTC)  headform, 
tested  in  different  cavitation  facilities  demonstrates  that 
cavitation  inception  on  a  given  model  can  have  many  dif¬ 
ferent  physical  forms  and  cavitation  inception  indices,  de¬ 
pending  on  the  environment  and  body  surface  condition.--'-' 
Departures  of  cavitation  inception  from  the  traditional  rule 
are  attributed  to  so-called  "scale  effects."  The  influence  of 
boundary  layer  pressure  fluctuations  on  cavitation  inception 
is  known  to  he  one  of  the  major  sources  of  scale  effect  -  1 


Rouse6  showed  that  the  high  pressure  fluctuations 
generated  in  the  turbulent  shear  layer  of  a  free-jet  can  pro¬ 
duce  cavitation  at  static  pressures  which  are  higher  than 
vapor  pressure.  Daily  and  Johnson8  showed  that  pressure 
fluctuations  in  the  middle  of  a  fully  established  turbulent 
boundary  layer  can  initiate  inception.  Levkovskii  and  Chalov7 
reported  that  turbulence  in  flow  past  a  wing  causes  cavitation 
to  incept  earlier.  Recent  work  by  Arakeri  and  Acosta"-9  and 
Huang4  indicates  that  for  supercritical  Reynolds  numbers  for 
which  bodies  do  not  have  laminar  boundary  layer  separation, 
cavitation  inception  takes  place  in  the  region  of  transition 
from  laminar  to  turbulent  flow.  For  subcritical  Reynolds 
numbers,  cavitation  inception  is  found  to  occur  in  the  reat¬ 
tachment  region  following  laminar  boundary  layer  separation. 

Significant  boundary  layer  pressure  fluctuations  on 
body  surfaces  have  been  measured  by  Huang  and  Hannan111 
and  by  Arakeri11  in  the  reattachment  region  following 
laminar  separation.  In  addition,  Huang  and  Hannan  reported 
that  measured  pressure  fluctuations  in  a  natural  transition 
region  are  higher  than  fluctuations  in  a  fully  established  tur¬ 
bulent  boundary  layer  flow  by  a  factor  ranging  from  2  to  3. 
Huang  and  Peterson12  reported  that  a  significant  scale  effect 
on  cavitation  inception  induced  by  boundary  layer  pressure 
fluctuations  exists  between  full-scale  and  model  propellers 
due  to  large  differences  in  Reynolds  numbers. 

The  properties  of  the  wall  fluctuating  pressure  field 
associated  with  the  growth  of  turbulent  spots  in  a  transition 
boundary  layer  on  a  flat  plate  have  been  measured  by  DeMetz 
et  al.,13  DeMetz  and  Casarella,'4  and  Gedney.15  Each  re¬ 
ported  that  pressure  fluctuations  in  the  transition  region  are 
smaller  than  those  in  the  fully  developed  turbulent  region. 
This  result  is  not  in  agreement  with  the  measurements  of 
Huang  and  Hannan.10 

A  basic  question  has  been  raised  as  to  whether  the 
boundary  layer  pressure  fluctuations  in  a  natural  transition 
region  can  be  greater  than  those  in  a  fully  established  tur¬ 
bulent  boundary  layer  region.  The  significance  of  this  point  is 
that  if  the  pressure  fluctuations  at  transition  are  greater  than 
in  the  fully  developed  turbulent  boundary  layer,  then  one  can 
expect  cavitation  to  preferentially  occur  at  the  transition  loca¬ 
tion.  In  view  of  the  significance  of  this  question  in  cavitation 
scaling,  fluctuating  pressure  measurements  fo:  a  two 
dimensional  hydrofoil  are  given  in  this  pape>- 

When  a  hydrofoil  operates  in  waves  or  a  propeller 
operates  behind  a  ship  hull,  the  lifting  surfaces  are  subjected 
to  temporally  and/or  spatially  varying  velocity  and  pressure 
fields.  The  time  varying  pressure  field  can  be  expected  to 
have  a  significant  effect  on  the  characteristics  of  cavitation 
inception,  growth,  and  collapse.  As  a  first  approximation,  the 
unsteady  effect  on  cavitation  has  been  studied  with  oscilla¬ 
ting  hydrofoils  by  Miyata  et  al.,lfi  Radhi,17  Bark  and  van 
Berlekom,1*  Shen  and  Peterson,19  and  van  Houten.20  The 
hydrofoil  approach  is  especially  attractive  for  simulating  a 
propeller  blade  operating  behind  an  inclined  shaft  with  the 
effective  angle  of  attack  oscillating  periodically  during  each 
cycle  of  rotation.  Available  data  demonstrate  that  the  cavi¬ 
tation-inception  angle-of-attack  at  the  leading  edge  is  differ¬ 
ent  from  inception  angles  determined  from  quasi-steady 
analysis  based  on  uniform  flow  experiments. 

Shen  and  Peterson19  showed  that  the  computed 
potential-flow  pressure  distributions  are  subject  to  a  signifi 
cant  phase  shift  with  respect  to  the  foil  oscillation  angle, 
which  in  turn  influences  the  occurrence  of  cavitation  incep¬ 
tion.  As  previously  mentioned,  significant  physical  effects  of 
boundary  layer  characteristics  on  cavitation  have  been  well 
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documented.  So  far  most  transition  and  boundary  layer  pres¬ 
sure  fluctuation  studies  have  been  limited  to  flat  plates  and 
headforms.  In  the  first  part  of  this  paper  the  influence  of 
oscillation  on  transition  boundary-layer  pressure  fluctuations, 
and  the  possible  effect  of  oscillation  on  cavitation  inception 
for  a  two-dimensional  hydrofoil,  will  be  presented.  As  pointed 
out  previously,  cavitation  frequently  is  unavoidable  on  the  lift¬ 
ing  surfaces  of  hydrofoils  and  marine  propellers.  In  the  sec¬ 
ond  part  of  this  paper  the  influence  of  oscillation  on  foil 
cavitation  instability  and  noise  will  be  presented. 


2.2  Test  Section 

The  closed-jet  test  section  of  the  36-inch  water  tunnel 
was  modified  by  the  insertion  of  sidewall  liners  to  provide 
two  flat  sides.  On  each  end  of  the  foil  a  disc  was  attached. 
This  disc  rotated  in  a  sidewall  recess;  thus  the  foil  could  be 
rotated  without  gap  cavitation  occurring  between  the  end  of 
the  foil  and  the  sidewall  of  the  tunnel.  One  sidewall  assembly 
was  fitted  with  clear  plastic  windows  to  permit  side  view 
photography. 


2.  EXPERIMENTAL  APPARATUS  AND  TEST 
PROCEDURE 

2.1  Foil  and  Instrumentation 


A  foil  was  machined  from  17-4  PH  stainless  steel  in 
the  form  of  a  rectangular  wing  of  Joukowski  section  with  the 
trailing  edge  modified  to  eliminate  the  cusp.  The  offsets  are 
given  by  the  following  equation: 


y/c 

e 

y/c 


0.04077  sin  0  -  0.02039  sin  (20) 
cos1  (1.888  x/c  -  1) 

0.08590  (1  -  x/c) 


0<  x/c  <0.7945 

(1) 

0.7945<  x/c<  1.0 


To  simulate  prototype  viscous  effects  at  the  leading  edge  as 
closely  as  possible,  die  model  was  designed  with  a  chord 
length  c  of  24.1  cm,  a  span  of  77.5  cm  and  a  relatively  large 
leading-edge  radius.  The  maximum  thickness  to  chord  ratio  is 
10.5  percent.  The  foil  surface  was  hand  finished  within 
0.38  pm  rms  surface  smoothness. 

Pressure  transducers  were  installed  at  a  distance  of 
7.96,  24.1  and  60.3  mm  from  the  leading  edge.  These  loca¬ 
tions  correspond  to  3.3,  10,  and  25  percent  of  chord  length 
from  the  leading  edge.  Kulite  semiconductor  pressure  gages 
of  the  diaphram  type  (Model  number  LQM- 10-250-305)  were 
mounted  within  a  Helmholtz  chamber  connected  to  the  foil 
surface  by  a  pinhole.  This  arrangement  permitted  measure¬ 
ment  of  unsteady  surface  pressures  due  to  foil  oscillation  and 
high  frequency  pressure  fluctuations  inside  the  boundary 
layer  over  a  pressure  range  of  ±  207  KPa  ( ±  30  PSI).  In  order 
to  increase  the  spatial  resolution  in  measuring  the  local 
pressure  fluctuations  inside  the  boundary  layer,  the  diameter 
of  the  pinholes  installed  on  the  foil  surface  was  kept  at  0.31 
mm  (0.012  inches).  This  arrangement  also  reduces  the  danger 
of  cavitation  damage  to  the  pressure  transducers.  Extreme 
care  was  taken  to  fill  the  Helmholtz-type  chamber  through 
the  pinhole  under  vacuum  with  deaerated  water  to  minimize 
the  possible  occurrence  of  an  air  bubble  trapped  inside  the 
chamber.  A  flat  pressure  response  was  observed  up  to  2,000 
Hz  in  dynamic  calibration  tests.  The  calibration  procedure 
used  here  was  developed  by  the  National  Bureau  of  Stand¬ 
ards,21  modified  to  the  extent  that  water,  rattier  than  silicone 
oil,  was  the  fluid  medium.  Since  it  was  very  important  to 
determine  the  relative  phase  difference  between  the  foil  angle 
and  the  pressure  gage  signals,  all  amplification  and  recording 
equipment  was  selected  to  minimize  the  introduction  of  un¬ 
wanted  phase  shifts. 

The  cavitation  noise  was  sensed  by  a  flush  mounted 
hydrophone  117  cm  downstream  from  the  foil  axis,  on  the  top 
of  the  water  tunnel's  close-jet  test  section.  The  hydrophone 
was  uncalibrated,  so  all  noise  measurements  are  relative  to  an 
arbitrary  level. 


Figure  1  —  Conceptual  Design  of  Foil  Oscillation  Mechanism 

The  foil  was  oscillated  in  pitching  motion  around  the 
quarter  chord  by  a  mechanism  whose  conceptual  design  is 
shown  in  Figure  1.  With  this  type  of  design,  the  foil  mean 
angle  (»o)  can  be  adjusted  statically  and  the  amplitude  of  foil 
oscillation  (oj)  can  be  continuously  adjusted  between  0°  to  4° 
while  in  operation.  The  oscillation  frequencies  cover  the  range 
between  4  Hz  to  25  Hz.  Air  bags  were  installed  to  reduce  the 
fluctuating  torque  requirements  on  the  motor  drive  system. 


2.3  Data  Reduction 

Due  to  installation  of  two  sidewall  liners  in  the  test 
section,  the  tunnel  speed  was  corrected  according  to  the  area- 
ratio  rule.  The  tape-recorded  time  histories  of  foil  angle  and 
pressures  were  digitized  by  using  a  Hewlett  Packard  2100 
minicomputer  and  reduced  by  using  algorithms  implemented 
on  the  DTNSRDC  CDC-6000  digital  computer.  The  time 
histories  were  recorded  on  one  inch  magnetic  tape  at  15  in 
ches  per  second  (38  cm/s)  using  IRK.'  standard  intermediate 
band  frequency  modulation  techniques,  (luring  digitization, 
pressure  fluctuation  data  were  filtered  using  four-pole  Butter- 
worth  bandpass  filters  that  have  a  -  3  dB  signal  attentuation 
at  80  and  2,000  Hz  for  Ixiundary  layer  pressure  fluctuations. 
The  run  lengths  used  in  the  data  reduction  were  nominally  40 
seconds.  For  the  oscillating  foil  data,  the  computer  output 
consisted  of  values  of  mean  and  standard  deviations  of  pres¬ 
sure  fluctuations. 
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3.  BOUNDARY  LAYER  PRESSURE  FLUCTUATIONS 
ON  A  STATIONARY  FOIL 

Consider  a  steady  uniform  flow  past  a  two-dimensional 
hydrofoil.  Let  the  local  pressure  p  on  the  foil  surface  be  ex¬ 
pressed  in  the  following  way: 

p  =  p,,  +  Ap,  (21 

where  p,,  is  the  mean  pressure  and  Ap,  is  the  pressure  assort- 
ated  with  the  lx>undary  layer  pressure  fluctuations.  Let 
denote  local  pressure  coefficient: 


where 


and 


P  -  Poo 

,;2  e  Voo* 

=  C ,,,  +  AC,, 


Pj.  ~  Poo 

“'2  e \V 


AC,„ 
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here  e,  Pco,  and  V*  are  fluid  density,  reference  free-stream 
static  pressure  and  reference  free-stream  velocity,  respec¬ 
tively. 

3.1  Laminar  Boundary  Layer  Stability  Calculations 


Figure  2  —  Sequence  of  Ensemble  Averaging 
(From  an  unpublished  report  by  R  Pierce,  DTNSRDC) 


The  hydrophone  signal  was  first  bandpass  filtered  bet¬ 
ween  10  kHz  and  40  kHz  to  minimize  signal  contamination 
from  noncavitation  related  sources.  This  signal  was  then 
squared  and  passed  through  an  averaging  filter  to  obtain  a 
“running  average,”  i.e.  quasi-stationary  sound  power  levels, 

p2  as  shown  in  Figure  2.  This  sound  power  level  is  also  digi¬ 
tized  along  with  the  signals  corresponding  to  the  foil  angle-of- 
attack  and  the  mean  foil  angle.  In  order  to  obtain  the  average 
relative  sound  power,  the  digitized  signal  is  averaged  over  the 
whole  run.  In  order  to  determine  the  average  waveform  of 
the  hydrophone  signal  during  one  cycle,  a  cycle  is  defined  as 
being  initiated  and  terminated  when  the  increasing  foil  angle 
passes  through  its  mean  position.  Each  cycle  is  divided  into 
(1600/f)  time  increments,  i.e.  at  a  f  =  4  Hz  oscillation  fre¬ 
quency,  400  time  increments  are  used.  The  digitized  hydro¬ 
phone  signal  in  each  time  increment  is  averaged  over  all  of 
the  cycles  in  a  run.  The  averages  in  each  of  the  increments 
then  collectively  represent  the  average  waveform  occurring 
for  one  cycle. 


Let  k  =  — —  denote  the  reduced  frequency  where  w  is 

^  V  oo 


the  circular  frequency  associated  with  foil  oscillation.  The 
limiting  case  k  =  0  corresponds  to  a  stationary  condition.  A 
uniform  flow  past  a  stationary  foil  at  an  angle-of-attack  a  of 
3.25  degrees  will  be  considered  first.  The  computed  steady 
potential  flow  pressure  distribution  for  this  foil  angle  is 
shown  in  Figure  3.  The  flow  field  to  be  measured  by  the 
pressure  gages  at  3.3,  10  and  25  percent  chord  length  is  seen 
to  be  in  a  strong  adverse  pressure  gradient.  This  fact  will  be 
shown  to  have  a  significant  effect  on  the  magnitude  of  bound¬ 
ary  layer  pressure  fluctuations.  The  prediction  of  laminar 
boundiary  layer  separation  is  based  on  the  criterion  of  zero 
shear  stress.  The  Smith  disturbance  amplification  method22  is 
used  to  correlate  and  predict  flow  transition. 

Let  A/Ao  denote  the  Smith  cumulative  spatial  amplifi¬ 
cation  ratio  A  to  its  amplitude  A,,  at  the  point  of  neutral  sta¬ 
bility.  Let  Rno  denote  the  Reynolds  number  based  on  chord 
length.  As  a  numerical  example,  the  computed  amplification 
ratio  versus  chordwise  location  xJc  is  given  in  Figure  4a  at 
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Figure  3  —  Measured  and  Predicted  Pressure  Distributions 

Rm  =  1.6  x  lCv  The  Smith  cumuiative  spatial  amplification 
ratio  was  computed  for  a  range  of  nondimensional  boundary 
layer  disturbance  frequencies  5, 


where  Q  is  the  dimensional  disturbance  frequency.  Figure  4a 
shows  that  disturbances  will  be  significantly  amplified  in  the 
frequency  range  of  Q  =  1.7  x  10 4  to  3.2  x  104  which  corres¬ 
pond  to  dimensional  frequencies  of  1,250  Hz_to  2,350  Hz. 
Figure  4a  also  indicates  that  the  frequency  O  =  3.2  x  10 4  has 
the  largest  amplification  ratio  up  to  x/c  =  0. 1 1 .  The  free- 
stream  velocity  used  in  these  computations  is  Va,  =  6.71  m/s. 
The  same  method  was  used  for  stability  calculations  at  other 
Reynolds  numbers.  The  computed  results  are  summarized  in 
Figure  4b. 

3.2  Experimental  Results 

The  experiments  with  the  foil  at  o  =  3.25  degrees 
were  performed  in  the  DTNSRDC  36-inch  water  tunnel.  As 
seen  in  Figure  3,  the  measurements  of  static  pressure  coeffi¬ 
cients  are  in  good  agreement  with  prediction.  At  speeds  of 
4.88,  6.71,  9.75,  13.11  and  14.94  meter/sec,  the  pressure 
responses  measured  at  three  foil  locations  are  given  in 
Figures  5a  to  5e.  The  pressure  gages  located  at  3.3,  10  and 
25  percent  chord  lengths  on  the  upper  surface  are  denoted  by 
P.t,  Pm  and  Pzs,  respectively.  The  purpose  of  steady  runs  was 


Figure  4b  —  Effect  of  Reynolds  Numbers  on  Amplification 
Ratios  at  a  =  3.25  deg 
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to  provide  reference  data  to  be  used  for  comparison  with 
dynamic  runs.  Consequently,  the  test  runs  given  in  Figures 
5a  to  5e  were  conducted  with  the  oscillating  mechanism  in 
motion  and  the  oscillating  amplitude  a i  set  to  zero  to  incor¬ 
porate  the  possible  effect  of  noise  and  vibration  on  the  meas¬ 
ured  phenomena. 

RUN  NO.  7001 


3  3  PERCENT  CHORD  LENGTH  p. 


10  PERCENT  CHORD  LENGTH  P10 


2S  PERCENT  CHORD  LENGTH  Pn 


Figure  5a  —  Measured  Pressure  Fluctuations  at  R„c  =  1.2  x  10* 
in  a  Steady  Flow  at  a  -  3.25 
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3  3  PERCENT  CHORD  LENGTH  P3 


10  PERCENT  CHORD  LENGTH  P,„ 

26  PERCENT  CHORD  LENGTH  PM 
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Figure  5d  —  Measured  Pressure  Fluctuations  at  Hnr  =  3.2  x  10,; 
in  a  Steady  Flow  at  a  =  3.25 
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3  3  PERCENT  CHORD  LENGTH  P3 


10  PERCENT  CHORD  LENGTH  P,0 


10  PERCENT  CHORD  LENGTH  P1( 


^  25  PERCENT  CHORD  LENGTH  Pn 

Figure  5b  —  Measured  Pressure  Fluctuations  at  R™  =  1.6  x  10* 
in  a  Steady  Flow  at  a  =  3.25 


RUN  NO.  8040 


3  3  PERCENT  CMO"D  LENGTH  P3 


10  PERCENT  CHORD  LENGTH  P10 

.  26  PERCENT  CHORD  LENGTH  Pa 

Figure  5c  —  Measured  Pressure  Fluctuations  at  RnC  -  2.4  x  10* 
in  a  Steady  Flow  at  a  »  3.25 


26  PERCENT  CHORD  LENGTH  Pa  t 

Figure  5e  —  Measured  Pressure  Fluctuations  at  R„c  =  3.6  x  10* 
in  a  Steady  Flow  at  a  =  3.25 

The  energy  spectrum  based  on  narrow  band  frequency 
analysis  from  80  to  5040  Hz  of  Run  8028  at  pressure  gage 
Pio  is  given  in  Figure  6.  The  maximum  fluctuating  signal  at 
Pio  corresponds  to  2650  Hz.  This  value  is  in  close  agreement 
with  the  computed  frequency  of  the  most  unstable  disturb¬ 
ances;  see  Figure  4a.  The  energy  spectra  of  pressure  responses 
based  on  narrow  band  frequency  analysis  were  conducted  for 
several  test  Runs.  Figure  6  shows  a  sharp  peak  at  4.7  kHz  at 
pressure  gage  Pio.  The  same  kind  of  sharp  peak  at  4.7  kHz 
was  observed  in  all  the  runs  analyzed.  A  similar  phenomenon 
was  observed  at  5.0  kHz  at  pressure  gage  P25.  It  is  suspected 
that  this  may  be  resonance  phenomenon  associated  with 
pinhole-type  pressure  gages.  The  pressure  gages  had  been 
dynamically  calibrated  up  to  2,000  Hz  with  flat  reponse.  Con¬ 
sequently,  the  data  given  in  Figures  5a  to  5e  were  bandpass 
filtered  with  4  pole  Butterworth  filters  from  80  to  2,000  Hz. 
Burton23  has  shown  that  the  energy  spectrum  of  a  flow  field 
in  an  adverse  pressure  gradient  is  compacted  in  a  much  nar¬ 
rower  and  lower  frequency  ranges  than  the  energy  spectrum 
associated  with  a  flow  field  in  a  zero  pressure  gradient. 

A  comparison  of  pressure  response  measurements  and 
stability  calculations  suggests  that  the  boundary  layer  at  P2;, 
(25  percent  chord  length)  for  all  the  runs  given  in  Figure  5a 
to  5e  is  fully  turbulent.  The  measured  root-mean-square  (rms) 
pressure  fluctuations  normalized  by  the  free-stream  dynamic 
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Figure  6  —  Energy  Spectrum  of  Run  8028  at  R^.  -  1.6  x  10* 

head  are  given  in  Table  1.  The  values  of  ACP,  range  between 
0.0065  to  0.011.  Burton  reported  that  the  measured  rms  pres¬ 
sure  fluctuation  coefficient  normalized  by  the  free-stream 
dynamic  head  is  measurably  smaller  in  an  adverse  pressure 
gradient  (around  0.0078)  than  in  a  favorable  pressure  gra¬ 
dient  (around  0.010).  Huang  and  Hannan10  measured  0.015  on 
an  axisymmetric  body.  The  present  data  were  digitized 
through  a  bandpass  filter  from  80  to  2,000  Hz.  The  loss  of 
high  frequency  range  may  influence  the  resolutions  of  pres¬ 
sure  gage  measurements  and  underestimate  the  magnitude  of 
pressure  fluctuations. 
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TABLE  1  -  ROOT  MEAN-SQUARE  BOUNDARY  LAYER 
PRESSURE  FLUCTUATIONS  AT  25  PERCENT  CHORD 
(P25)  AND  10  PERCENT  CHORD  (Pln) 


Additionally,  Table  1  shows  that  at  a  given  location 
the  measured  pressure  fluctuation  coefficients  decrease  with 
an  increase  in  speed.  The  same  trend  was  also  reported  by 
Burton.  Measurements  by  Ludwieg  and  Tillman24  showed 
that  the  shearing  stress  in  a  turbulent  boundary  layer  de¬ 
creased  significantly  with  an  increase  in  momentum  thickness 
Reynolds  number.  Blake25  reported  that  the  shear  stress 
rather  than  the  free-stream  dynamic  head  should  be  used  to 
collapse  the  data.  All  these  findings  suggest  that  the  bound¬ 
ary-layer  pressure  fluctuation  coefficient  to  be  used  for 
model/full-Bcale  cavitation  scaling  is  not  a  universal  constant, 
as  used  in  Reference  12,  but  instead  is  a  function  of  Reynolds 
number. 


Stability  calculations  at  Rn<.  -  1.6  x  106  show  that 
amplification  ratios  of  e®,  es  and  e11  are  reached  at  xJc  - 
0.10,  0.13  and  0.16,  respectively.  Laminar  separation  is 
predicted  to  occur  at  x/c  -  0.23.  Theoretical  computations 
suggest  that  natural  transition  from  laminar  to  turbulent 
boundary  layers  is  to  be  expected  to  occur  upstream  of  the 
laminar  separation  point  at  this  Reynolds  number.  Thus, 
natural  transition  precludes  the  possibility  of  laminar  separa¬ 
tion. 

At  a  speed  of  4.88  m/s,  Figure  5a  shows  that  the  flow 
at  3.3  and  10  percent  chord  lengths  are  laminar.  With  an  in¬ 
crease  in  speed  to  6.71  m/s,  corresponding  to  R^  =  1.6  x  10®, 
a  strong  pressure  fluctuating  signal  with  intermittency  is 
noticed  at  10  percent  chord;  see  Figure  5b.  The  bursting 
signal  observed  in  Figure  5b  is  related  to  natural  transition. 

It  was  found  that  the  computed  amplification  ratio  of  e®  cor¬ 
related  well  with  the  measured  location  of  natural  transition 
for  the  present  experimental  setup  at  Rn,  =  1.6  x  10®.  For 
other  Reynolds  numbers  the  computed  locations  of  the 
amplification  ratios  of  e®  were  used  to  indicate  the  start  of 
natural  transition 

At  a  speed  of  9.75  m/s.  corresponding  to  R,,.  =  2.4  x 
10®,  the  stability  calculations  predict  an  amplification  ratio  of 
e85  at  P ) , Measured  pressure  responses  suggest  that  the 
flow  at  10  percent  chord  is  in  the  final  phase  of  transition 
with  an  intermittency  factor  close  to  unity;  see  Figure  5c. 

The  rms  pressure  fluctuations  normalized  by  the  free-stream 
dynamic  head,  are  given  in  Table  1.  The  measured  rms  pres¬ 
sure  fluctuation  at  10  percent  chord  in  the  transition  region 
are  seen  to  be  greater  than  those  measured  at  25  percent 
chord  in  the  turbulent  region.  This  observation  is  based  on 
data  analysis  with  a  bandpass  filter  of  80  to  2,oou  Hz.  In  the 
few  cases  when  the  bandpass  filter  upper  limit  increased  to 
5,040  Hz,  the  same  conclusion  held. 

The  properties  of  the  wall  fluctuating  pressure  field 
associated  with  the  growth  of  turbulent  spots  in  a  natural 
transition  boundary  layer  on  a  flat  plate  have  l>een  measured 
by  DeMetz  et  al.,l:l  DeMetz  and  Uasarella,14  and  Uednev.1 ' 
They  all  report  that  the  pressure  amplitudes  in  the  wall 
pressure  field  of  individual  spots  are  approximately  equal  to 
those  in  the  fully  turbulent  boundary  layer  region.  As  the  in¬ 
termittency  factor  y  increases,  the  power  sjiectral  densities 
increase  systematically  to  a  maximum  value  in  the  fully  tur¬ 
bulent  condition.1'' 

p-'(y)  =  rp-'(r),  ,  ;  (7| 

These  results  suggest  that  the  pressure  fluctuations  measured 
in  a  natural  transition  region  are  smaller  than  those  in  a  tur¬ 
bulent  region,  which  is  not  in  agreement  with  the  present 
measurements. 

Pressure  fluctuations  in  regions  of  natural  flow  transi¬ 
tion  on  an  axisymmetncal  body  were  also  measured  by  Huang 
and  Hannan.1"  Compared  to  the  pressure  flucturations  meas¬ 
ured  in  a  fully  established  turbulent  region,  the  intensity  of 
pressure  fluctuations  in  the  transition  region  is  higher  thar  in 
the  turbulent  region.  The  present  results  agree  with  the 
results  of  Huang  and  Hannan  These  measurements  and  the 
present  measurements  were  made  in  flows  with  adverse  pres 
sure  gradients,  whereas  the  linear  relationship  given  in  Equa¬ 
tion  (7)  is  derived  from  measurements  on  a  flat  plate  with 
zero  pressure  gradient.  It  is  suspected  that  the  cause  of  dif¬ 
ferences  in  the  two  sets  of  data  is  differences  in  the  pressure 
gradients.  This  difference  is  further  amplified  in  the  following 
series  of  tests  The  pressure  fluctuations  at  a  speed  of  13  11 
m/s,  corresponding  to  R,„.  =  3.2  x  10®.  are  given  in  Figure 


4.2  Flow  Description 


5d.  This  figure  suggests  that  the  flow  is  turbulent  at  10  per¬ 
cent  chord;  this  correlates  with  the  stability  calculations 
which  indicated  that  the  amplification  ratio  reached  about  e9 
to  e10  at  P,0.  Figure  Sd  also  suggests  that  the  location  of 
transition  has  moved  ahead  of  Pio-  The  measured  pressure 
coefficient  at  this  speed  is  smaller  than  the  one  measured  at  a 
speed  of  9.75  m/s.  This  fact  seems  to  support  the  argument 
that  the  pressure  fluctuation  in  a  transition  region  is  greater 

than  in  a  turbulent  region  with  an  adverse  pressure  gradient. 
A  further  increase  in  speed  to  14.94  m/s  reduces  the  pressure 
fluctuation  further  as  shown  in  Figure  5e. 

As  a  final  remark.  Figure  5d  shows  that  the  frequent¬ 
ly-occurring  large  negative  pressure  fluctuations  are  approx¬ 
imately  2.5  times  higher  than  their  rms  value  in  the  transition 
region  at  P1(,  and  2.0  times  higher  than  their  rms  value  in  the 
fully  established  turbulent  region  at  P2s.  These  results  are  in 
close  agreement  with  the  Huang  and  Hannan  measure¬ 
ments.!0  Huang  and  Peterson12  suggested  that  the  frequent¬ 
ly-occurring  large  negative  pressure  fluctuations  in  transition 
may  be  responsible  for  cavitation  inception  events. 


4.  BOUNDARY  LAYER  PRESSURE  FLUCTUATIONS 
WITH  FOIL  IN  OSCILLATION 

4.1  The  Oscillating  Foil 

The  foil  oscillates  sinusoidally  in  pitching  motion 
around  the  quarter  chord  point,  measured  from  the  foil 
leading  edge.  The  instantaneous  foil  angle  is  denoted  by: 


For  a  typical  propeller  operating  behind  an  inclined 
shaft,  the  lifting  surface  representing  a  propeller  blade  will 
encounter  a  periodic  pressure  field.  The  amplitude  of  pressure 
fluctuation  depends  on  the  ship  speed  and  the  angle  of  in¬ 
clination.  Angle-of-attack  variations  of  0.5  to  2.0  degrees  at 
the  0.7  propeller  radius  are  expected  to  represent  a  range  of 
practical  interest.  The  effect  of  unsteady  angle-of-attack 
variations  on  lifting  surface  boundary-layer  pressure  fluctua¬ 
tions  was  investigated  by  pitching  the  previously -described 
two-dimensional  hydrofoil  sinusoidally. 

A  typical  oscillograph  plot  of  pressure  responses  with 
the  foil  in  oscillation  is  shown  in  Figure  7.  The  data  were 
bandpass  filtered  from  80  to  2,000  Hz.  The  low  frequency 
periodic  part  has  been  subtracted.  To  accommodate  an  oscilla¬ 
tion  cycle  of  the  foil  on  a  reasonably  sized  sheet,  the  trace 
was  processed  at  2  inches  per  second  paper  speed.  (Figures 
5a  to  5e  were  processed  at  10  inches  per  second  paper  speed.) 
The  first  line  from  the  top  gives  the  instantaneous  foil  angle. 

In  this  example,  the  frequency  of  foil  oscillation  is  4  Hz.  The 
instantaneous  foil  angle  is  given  by  o  =  3.25°  +  2.  ] 0  sin  wt. 
The  speed  is  9.75  m/s,  corresponding  to  R„r  *  2.4  x  10*. 

Recall  that  in  the  steady  run  (a  »  3.25°)  at  R,,,  =  2.4  x  10* 
(Figure  5b),  the  flow  was  in  the  final  stage  of  transition  at 
the  10  percent  chord  length.  The  low  amplitude  unsteady 
signal  appearing  at  3.3  percent  chord  length  is  due  to 
mechanical  noise  caused  by  the  oscillator  motion  and  the  elec¬ 
trical  noise  within  the  recording  electronics. 


o  =  o„  +  o,  sin  (cot) 


(8) 
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where  <r0,  a\,  and  co  are  the  mean  foil  angle,  amplitude  of  fod 
oscillation  and  oscillation  frequency,  respectively.  Let  ♦ 
denote  the  phase  angle  between  the  pressure  response  and 
the  instantaneous  foil  angle.  At  a  given  location  on  the  foil, 
the  time  dependent  local  pressure  and  pressure  coefficient 
are  expressed  by: 


P(t)  >  P,  +  |  Pu  sin  (cot  +  4)  +  AP, 

Cp(t)  **  Cps  +  |  ACpU|  sin  (cot  +  4)  +  ACP, 

where 


Ps  -  Poo 

1/2  C  V 


1 4CpU| 


I  Pul 

V*  e  V*2 
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Figure  7  -  Sample  Oscillograph  Record  at  R,„  =  2  4  x  10\ 
f  =  4  Hz  and  o  »  3.25  +  2. 1  sin  cot 


Here,  Cp,,  |  ACpu|  and  ACP,  respectively  denote  mean  poten¬ 
tial  flow  pressure  coefficient,  dynamic  potential  flow  pressure 
coefficient,  and  the  unsteady  pressure  coefficient  associated 
with  boundary  layer  pressure  fluctuations.  The  objective  of 
the  work  described  here  was  to  determine  the  influence  of  foil 
oscillation  on  transition  and  magnitudes  of  boundary  layer 
pressure  fluctuations. 


In  order  to  describe  the  observed  chahge  in  transition 
location  with  foil  angle,  consider  the  pressure  response  pt  25 
percent  chord  length,  P-:„  Figure  7  suggests  that  the  flow  is 
turbulent  and  that  the  location  of  transition  is  ahead  of  p . 
when  the  foil  is  at  the  maximum  foil  angle  a  =  onmx  As  the 
foil  angle  begins  to  decrease,  the  location  of  transition  moves 
toward  the  trailing  edge.  At  some  intermediate  foil  angle  o  = 
»i,  the  signals  of  pressure  fluctuations  are  intensified  and  the 


transition  is  assumed  to  occur  at  P^.  After  the  foil  reaches  a 
-  the  foil  angle  begins  to  increase  and  the  movement  of 
transition  is  reversed  toward  the  leading  edge.  Finally,  the 
flow  at  becomes  fully  turbulent  again  at  a  -  or- 

Next,  consider  the  pressure  response  at  10  percent 
chord  length,  Pio.  In  this  example,  the  minimum  foil  angle  is 
»mui  »  1.15°.  Figure  7  suggests  that  the  flow  is  laminar  at 
Pio  when  the  foil  angle  is  at  om^.  The  location  of  transition  is 
aft  of  the  10  percent  chord  length.  As  the  foil  angle  begins  to 
increase,  transition  moves  forward  with  instability  occurring 
at  a  certain  foil  angle,  a  »  <n.  A  further  increase  in  the  foil 
angle  produces  full  turbulence  at  a  »  ay.  It  should  be 
remarked  that  the  pressures  at  Pio  must  be  multiplied  by  a 
factor  of  1.15  to  account  for  difference  in  gage  sensitivities  if 
a  direct  comparison  of  pressure  fluctuations  between  P25  and 
Pio  are  to  be  made  on  the  oscillograph  plot.  The  flow  remains 
turbulent  when  the  foil  angle  reaches  o  -  This  is  follow¬ 
ed  by  a  decrease  in  foil  angle.  At  a  «  ay,  flow  transition  oc¬ 
curs  again  at  Pio.  A  further  decrease  in  foil  angle  produces 
laminar  flow  at  a  «  ol-  The  rms  pressures  at  Pio  derived 
from  Figure  7  are  given  in  Figure  8.  Figure  8  shows  that  the 
rms  value  of  4Cpt  attains  the  maximum  of  0.0155  at  the  foil 
angle  a  -  2.6  degrees.  Figure  7  gives  a  measured  large 
negative  pressure  fluctuations  of  AC,„  »  0.042,  which  is  ap¬ 
proximately  2.7  times  higher  than  the  peak  value  (rms)  shown 
in  Figure  8.  This  result  is  in  agreement  with  the  previous  sta¬ 
tionary  runs. 
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Recall  the  case  o>  «  0,  corresponding  to  Run  8040  as 
given  in  Figure  5c.  The  flow  at  P3  (3.3  percent  chord  length) 
is  seen  to  be  laminar.  The  flow  at  Pn,  is  seen  to  be  in  the 
final  phase  of  transition.  The  flow  at  Pzs  is  seen  to  be  fully 
turbulent. 

4.3.1  Pressure  Fluctuations  at  10  Percent 

Chord  Length  (Pm) 

The  case  with  oi  =  0.5  degrees  is  considered  first.  It 
is  remarked  that  the  flow  at  3.3  percent  chord  is  laminar 
throughout  the  whole  series  of  runs.  Because  the  pressure 
fluctuations  in  dynamic  runs  are  non-stationary,  the  test  data 
were  digitized  to  obtain  a  rms  value  (see  Figure  2).  The  cor¬ 
responding  rms  pressure  fluctuations  at  P)(,  are  shown  in 
Figure  9  for  oscillation  frequencies  of  4,  5.5  and  15  Hz.  The 
vertical  axis  gives  rms  pressure  fluctuations  normalized  by 
the  free-stream  dynamic  head.  The  horizontal  axis  gives  the 
instantaneous  foil  angle  a  and  the  oscillation  angle  wt.  The 
magnitudes  of  the  pressure  fluctuations  are  seen  to  vary  with 
the  instantaneous  foil  angles.  Figures  5c  and  9  suggest  that 
the  location  of  transition  is  slightly  aft  of  P,,,  at  a  =  3.25 
degrees.  As  the  foil  angles  are  increased,  transition  moves 
forward  toward  Pm.  The  measured  rms  pressure  fluctuations 
at  Pm  are  systematically  intensified.  As  the  foil  angle  ap¬ 
proaches  3.60  degrees,  the  location  of  transition  occurs  at  Pjo 
and  the  magnitude  of  pressure  fluctuations  attains  a  max¬ 
imum.  A  further  increase  in  foil  angle  moves  the  transition 
ahead  of  Pio;  and  the  magnitude  of  pressure  fluctuations 
measured  at  Pj0  is  reduced.  After  the  foil  angle  reaches  om„ 
the  foil  angle  begins  to  decrease,  transition  moves  down¬ 
stream  toward  Pio  and  the  intensity  of  pressure  fluctuations 
measured  at  P]0  increases  again.  At  a  foil  angle  of  about  3.25 
degrees  (wt  «  180  degrees),  the  location  of  transition  approx¬ 
imately  coincides  with  Pi0.  A  further  decrease  in  foil  angle 
moves  the  transition  location  toward  the  trailing  edge  and  the 
pressure  fluctuations  measured  at  P1M  begin  to  decrease. 
Transition  is  seen  to  pass  through  Pm  twice  in  every’  cycle  of 
oscillation.  This  results  in  the  appearance  of  two  peaks  in 
pressure  fluctuations.  The  whole  sequence  of  variation  in 
pressure  fluctuations  with  instantaneous  foil  angle  can  l>e 
best  illustrated  with  large  oscillation  amplitudes. 

*  > 


Figure  8  —  Measured  RMS  Pressure  Fluctuations  on  Run  8043 
at  Rnr  =  2.4  x  106  and  a  »  3.25  +  2.1  sin  wt 


4.3  Unsteady  Effects 

A  series  of  dynamic  runs  was  conducted  to  examine 
unsteady  effects  on  boundary  layer  pressure  fluctuations.  A 
free-stream  speed  of  V*,  *  9.75  m/s  was  used  for  this  series 
of  experiments.  The  instantaneous  foil  angle  o  is  given  by  a 
=  3.25  ♦  oi  sin  wt.  The  oscillation  amplitude  a\  covered  the 
range  from  0  to  2.1  degrees.  The  oscillation  frequency,  f  = 
w/2n,  covered  the  range  from  4  to  15  Hz.  Due  to  the  in¬ 
creased  background  noise,  the  pressure  fluctuation  associated 
with  the  25  Hz  foil  oscillation  are  not  presented  here. 
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Figure  9  —  Measured  RMS  Pressure  Fluctuations  at  l’j„  with 
Oscillating  Amplitude  of  0  5  Degrees  and  R,„  =  2  4  x  lu” 
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Figure  10  —  Measured  RMS  Pressure  Fluctuations  at  Pi0  with 
Oscillating  Amplitude  of  0.97  Degrees  and  R^  -  2.4  x  106 
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Figure  11  —  Measured  RMS  Pressure  Fluctuations  at  Pu>  with 
Oscillating  Amplitude  of  1.57  Degrees  and  -  2  4  x  106 
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Figure  12  —  Measured  RMS  Pressure  Fluctuations  at  Phi  with 
Oscillating  Amplitude  of  2  1  Degrees  and  R,„  -  2.4  x  10* 


1  At  the  10-percent  point  on  the  foil,  there  are  two 
peaks  of  pressure  fluctuations  associated  w  .th  ad 
vancing  and  receding  transition  locations  during 
each  cycle  of  oscillation  The  magnitudes  of  the  rms 
pressure  amplitudes  for  both  advancing  and  reced 
ing  boundary  layer  transitions  are  the  sanii-  at  low 
values  of  oscillation  frequency  14  and  5  f>  Hz) 
However,  at  high  values  of  oscillation  frequency  (111 
and  15  Hz)  the  magnitude  associated  with  ad vane 
ing  transition  is  generally  greater  than  the  magni 
tude  of  the  receding  transition 

2  The  peak  values  of  rms  pressure  fluctuations 
associated  with  advancing  transition  are  seen  to  lie 
independent  of  oscillation  amplitudes  and  frequen 
cies  and  they  are  almost  the  same  with  the  values 
measured  under  steady  conditions  These  results 
provide  a  contrast  between  the  Umndary  layer 
pressure  fluctuations  and  unsteady  potential  flow 
pressure  distributions  In  the  previous  )«i[ier  '  ‘  »e 
showed  that  the  magnitudes  of  dynamic  pressure 
responses  of  the  p4>tential  flow  (first  haTtinfnci  im¬ 
pend  linearly  on  the  oscillation  amplitudes  and  vary 
significantly  with  oscillation  frequencies 

3.  In  the  previous  paper,1*  we  were  able  to  successful 
ly  correlate  all  of  the  data  on  the  phase  relationship 
between  the  peak  of  dynamic  pressure  response  in 
the  time-varying  mean  flow  and  the  instantaneous 
foil  angle  by  use  of  the  reduced  frequency  param 
eter  K.  However,  in  the  present  pa[*-r.  due  to  the 
scattering  of  the  data,  we  can  only  qualitatively 
state  that  the  locations  of  boundary  layer  transition 
are  delayed  with  an  increase  in  reduced  frequency 

4  .  Consider  pressure  fluctuations  in  the  region  t>et 

ween  90<  w<  180  degrees  It  wu  stated  previous 
ly  that  in  this  region,  transition  is  located  forward 
of  P],i  and  that  the  magnitude  of  pressure  fluctua 
tions  decreases  systematically  with  increase  of 
oscillation  amplitude  Referring  to  Figure  12  for  an 
oscillation  amplitude  »,  >  2.10  degrees  the  fart 
that  the  pressure  fluctuations  remain  at  a  constant 
value  of  ACp,  *  0.005  for  90  <  uit  <  180,  indicates 
that  the  flow  is  fully  turbulent  in  this  region  This 
value  (ACp,  a  0.005)  measured  at  P,..  is  smaller 
than  the  value  (AC,„  a  0  008)  measured  at  P“,  with 
turbulent  flows  at  both  locations  It  is  not  dear 
whether  the  magnitudes  of  adverse  pressure  gra 
dients  in  these  two  locations  control  that  difference 

5  .  Consider  the  flow  characteristics  in  the  region 

around  270 <  uit  <  300  degrees  It  was  stati-d 
previously  that  the  transition  [mint  is  located  fai  aft 
of  P,„.  With  an  increase  in  amplitude  "f  oscillation 
the  flow  becomes  laminar  again  at  P  ,  as  seen  m 
Figures  11  ami  12  where  the  AC,,,  signal  observed 
in  the  range  27n<  uit<  lUKi  may  lie  attributeii  h> 
noise  contamination 


Figures  10,  11,  and  12  give  the  rms  pressure  fluctua¬ 
tions  measured  at  Pj,,  for  pitch  oscillation  amplitudes  of  0.97, 
1  57  and  2.10  degrees.  The  general  patterns  among  these 
figures  are  similar  Based  on  this  series  of  experiments,  we 
can  make  the  following  observations  concerning  pressure  flue 
tuations  at  the  10  percent  chord  location: 


1  > 


2  ff 

Iff 


fOU  ANGtl 
)  »  '<*•«! 


Figure  13  —  Measured  RMS  Pressure  Fluctuations  at  P .  with 
Oscillating  Amplitude  of  0.97  Degrees  and  R,„  *  2  4  x  lie 
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Figure  14  -  Measured  RMS  Pressure  Fluctuations  at  Po,  with 
Oscillating  Amplitude  of  1  57  Degrees  and  R,,  •-  2.4  x  10* 
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Figure  15  -  Measured  RMS  Pressure  Fluctuations  at  P25  with 
Oscillating  Amplitude  of  2.1  Degrees  and  R*.  •  2  4  x  10* 

4  .'1  2  Pressure  Fluctuations  at  25  Perce-t  Chord  length  (Pe.) 

Figures  13,  14,  and  15  give  rms  pressure  fluctuations 
at  PJ&  for  oscillating  amplitudes  of  0.97,  1.57,  and  2  10 
degrees,  respectively  The  flow  at  P2f,  is  fully  turbulent  with 
oi  •  0  5  degrees  (not  shown  in  the  figures)  throughout  each 
cycle  of  oarillation  Baaed  on  this  senes  of  experiments,  we 
can  make  the  following  observations 


1  The  value  of  rms  pressure  fluctuations  in  a  fully 
established  turbulent  flow  (ACpt  ■=  0  0080)  is  in¬ 
dependent  of  oscillating  amplitudes  and  frequencies. 
This  value  is  almost  identical  to  the  value  measured 
in  the  steady  runs 

2  With  an  increase  in  oscillation  amplitude,  transition 
moves  toward  P2j,  in  the  region  270  <U  4  360  The 
magnitudes  of  pressure  fluctuations  are  consequent¬ 
ly  intensified  at  P2«,  as  seen  in  Figure  15  Simultan¬ 
eously,  the  flow  becomes  laminar  at  P,,,  (see  Figure 
12!  Consistent  with  the  observation  at  P|(l,  the  oc¬ 
currence  of  peak  pressure  Auctions  at  P2r,  is 
delayed  by  increasing  the  oscillation  frequency 

3  The  influence  of  advancing  and  receding  transition 
points  on  pressure  fluctuations  measured  at  P2f,  is 
seen  to  be  compatible  with  the  behavior  at  P.,, 

4  ,'xs  shown  in  Figures  11  to  15  Ixiundary  layer  tran 
sition  in  an  oscillating  foil  is  delayed  by  an  increase 
in  oscillation  frequency  It  is  further  noticed  that 
the  pressure  fluctuations  in  a  transition  region 
(AC.,1  *  d  016)  are  greater  than  in  a  fully  turbulent 
region  (A('pt  *  0  OOHOI,  by  a  factor  ,,f  2  This  result 
is  in  agreement  with  the  observations  in  steady- 
runs 

4  4  Reynolds  Number  Kffects 

Consider  a  local  Reynolds  numlier  R,„  based  on  chord 
length  at  the  0  7  propeller  radius  Because  the  local  Reynolds 
number  associated  with  a  mode)  propeller  may  vary  from  one 
test  facility  to  another,  additional  dynamic  runs  were  con¬ 
ducted  to  examine  the  effect  of  Reynolds  number  on  bound 
ary  layer  pressure  fluctuations  These  additional  runs  were 
made  at  R„,  values  of  1  2  x  10*.  16x10*  and  3  6  x  16*  The 
earlier  runs,  a! read i  discussed,  were  made  at  the  R„,  values 
of  2.4  x  10* 

4  4  1  Pressure  Fluctuations  at  10- Percent  Chord  length  (P|  ,1 

Figures  16  and  17  give  the  rms  pressure  fluctuations 
measured  at  P,„  for  R„,  «  1  6  x  10*  anil  1  2  x  10*.  respective 
ly  The  examples  given  in  these  two  figures  are  for  an  oscilla 
tion  amplitude  of  a,  -  2  10"  Once  again,  two  i>eaks  in 
pressure  fluctuations  associated  with  advancing  and  receding 
transition  locations  are  observed  for  each  cycle  of  oscillation 
The  magnitudes  of  these  two  peaks  are  markedly  different 
Additionally,  the  rms  magnitudes  associated  with  advancing 
transition  are  significantly  intensified  at  the  lowest  value  of 
Rnr  -  1.2  x  10*.  reaching  as  high  as  10  percent  of  the 
dynamic  head  The  ordinate  in  Figure  1 7  has  been  modified  to 
accommodate  the  significant  increase  in  magnitude 

The  values  of  cumulative  spatial  amplification  ratio 
and  locations  of  laminar  separation  computed  for  steady  flow 
are  shown  in  Figure  18  It  was  previously  shown  in  steady 
runs  that  the  bursting  signal  associated  with  natural  transi 
tion  in  the  present  test  set-up  is  best  correlated  with  an 
amplification  ratio  of  e*  The  ohordwise  locations  for  di  >ti  r 
hances  to  achieve  amplification  ratios  of  e*  are  shown  in 
Figure  18  for  three  values  of  Reynolds  number  Within  the 
range  of  foil  angles  tested,  laminar  boundary  layer  stability- 
computations  indicate  that  natural  transition  will  occur  earlier 

than  laminar  separation  when  R,„  *  2  4  x  III*  On  the  other 
band,  laminar  separation  will  occur  prior  to  natural  transition 
at  foil  angles  aliove  4  H  degrees  at  R,  ,  1  6  x  111*  and  4  2 

degrees  at  K. ,  *  I  2  x  I  <  r 
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rure  16  —  Measured  RMS  Pressure  Fluctuations  at  Pjq  with 
Oscillating  Amplitude  of  2.1  Degrees  and  R^  -  1.6  x  10* 
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Figure  17  -  Measured  RMS  Pressure  Fluctuations  at  Pi0  with 
Oscillating  Amplitude  of  2.1  Degrees  and  R„  -  1.2  x  10* 


Figure  18  —  Chordwise  Location  of  Transition  and  Laminar 
Separation  Versus  a  and  R„(. 


In  comparison.  Figure  12  shows  the  values  of  AC  .,  at 
an  oscillation  amplitude  of  2.1  degrees  and  ICk  »  2.4  x  lo*. 
The  maximum  pressure  fluctuations  associated  with  advanc¬ 
ing  transition  occur  at  about  a  »  4.5  degrees  for  an  oscilla¬ 
tion  frequency  of  5.5  Hz.  Figures  4  and  18  suggest  that  the 
maximum  pressure  fluctuations  are  related  to  a  natural  tran 
sition  phenomenon.  The  rms  pressure  fluctuations  due  to  ad¬ 
vancing  transition  are  ACpi  -  0.016.  On  the  other  hand,  at 
the  lowest  value  of  R«  examained,  Rnc  -  1.2  x  10*.  Figure  17 
shows  that  the  maximum  pressure  fluctuations  occur  at  about 
5.3  degrees  at  5.5  Hz  with  rms  pressure  fluctuations  reaching 
to  ACp<  »  0.10.  Figure  18  suggests  that  this  peak  pressure  is 
associated  with  laminar  separation.  The  same  discussion  is  ap¬ 
plicable  to  Figure  16.  The  present  result  is  in  agreement  with 
Huang  and  Hannan's  finding  that  pressure  fluctuations 
associated  with  laminar  separation  are  much  greater  than  for 
natural  transition. 

Figures  12,  16  and  17  show  that  the  occurrence  of  ad 
vancing  transition  sets  in  earlier  as  Reynolds  number  in¬ 
creases.  This  trend  is  in  agreement  with  the  boundary  layer 
stability  calculations  shown  in  Figure  4  To  demonstrate  this 
fact  more  clearly,  the  pressure  fluctuations  measured  from 
three  dynamic  runs  at  R^  -  3.6  x  10*.  2  4x10*  and  1  2  x 
10*  are  shown  in  Figure  19  The  oscillation  amplitude  and  fre- 
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quency  associated  with  this  example  are  0.97  degrees  and  10 
Hz,  respectively.  For  -  1.2  x  106,  peak  pressure  fluctua¬ 
tions  never  reaches  Pio.  Consequently,  only  one  major  peak 
pressure  associated  with  laminar  separation  is  noticed  in  the 
whole  cycle  of  oscillation.  In  Figure  19  the  associated  back¬ 
ground  noise  was  corrected  in  such  a  way  that  ACpt  became 
zero  in  the  laminar  flow  region.  Due  to  the  phase  shift 
associated  with  the  potential  flow  pressure  distributions  and 
the  foil  angles,  the  process  of  becoming  laminar  flow  is 
achieved  at  a  foil  angle  greater  than 


Figure  19  —  Measured  RMS  Pressure  Fluctuations  at  P|»  with 
Three  Reynolds  Numbers  and  Oscillating  Amplitude  of 
0.97  Degrees  (f  -  10  Hz) 


4.4.2  Pressure  Fluctuations  at  25  Percent  Chord  Length  (P^) 

Figures  20  and  21  give  rms  pressure  fluctuations 
measured  at  P2S  for  R«r  »  16x  10*  and  1.2  x  10*.  respective¬ 
ly.  These  two  figures  support  the  previous  observation  that 
ACp,  in  a  fully  established  turbulent  region  is  independent  of 
oscillation  frequency.  The  same  conclusion  may  be  drawn 
from  Figure  22.  Additionally,  the  values  of  AC,,,  are  seen  to 
increase  with  decreasing  Reynolds  numbers. 
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Figure  20  —  Measured  RMS  Pressure  Fluctuations  at  P2s  with 
Oscillating  Amplitude  of  2.1  Degrees  and  -  1.6  x  10* 


Figure  21  —  Measured  RMS  Pressure  Fluctuations  at  P25  with 
Oscillating  Amplitude  of  2. 1  Degrees  and  R^  -  1.2  x  10* 
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Figure  22  —  Measured  RMS  Pressure  Fluctuations  at  P25  with 
Three  Reynolds  Numbers  and  Oscillating  Amplitude  of 
0.97  Degrees  (f  -  10  Hz) 


5.  INFLUENCE  OF  HYDROFOIL  OSCILLATION  ON 
CAVITATION  INCEPTION 


Let  o  denote  the  vapor  cavitation  number: 

P»  -  Pv 
°  '  WeV„2 


00) 


where  P,  is  the  vapor  pressure.  If  thermodynamic  equilibrium 
prevails,  then 


•  -  -Cp™,  (ID 

at  cavitation  inception  in  pure  fluids.2  In  practice  significant 
deviation  from  this  equation  has  been  widely  observed  in  the 
model  cavitation  testa.  Arndt2*  stated  that  the  so-called  scale 
effects  are  due  to  deviations  in  two  basic  assumptions  in¬ 
herent  in  the  cavitation  scaling  law,  namely  that  the  pressure 
scales  with  velocity  squared  and  the  critical  pressure  for  in¬ 
ception  is  the  vapor  pressure.  In  practice,  these  two  factors 
are  interrelated,  since  the  critical  pressure  is  a  function  of  the 
time  scale  of  the  pressure  field.  Consequently,  depending  on 
the  cavitation  resistance  of  the  liquid,  cavitation  inception  on 
the  model  may  occur  either  (a)  at  the  location  of  CPmin,  or  (b) 
in  the  natural  transition  or  laminar  separation  region  of  the 
model  12 

5  1  Cavitation  Inception  at  C^ 

We  will  first  consider  the  case  when  cavitation  incep¬ 
tion  on  the  model  occurs  at  the  location  of  Cpmin  Since  the 
magnitude  of  pressure  fluctuations  associated  with  a  laminar 
flow  is  esaentudly  zero,  the  present  experiments  show  that 


the  flow  at  3.3  percent  chord  length  is  laminar  for  ail  of  the 
dynamic  runs  so  that  ACpt  -  0.  Theoretical  computations  a  - 
3.25°  show  that  Cp^  occurs  around  x/c  -  0.018;  see  Figure 
3.  Thus,  Equation  (9)  reduces  to: 


Cp(t) 


Up*  +  Upu(t) 

Cp^  +  |  ACpul  sin  (wt  ♦ 


(12) 

for  0<  x/c  S  0.033 


where  |ACpu|  and  4  are  the  amplitude  and  phase  angle. 


An  earlier  series  of  experiments  conducted  by  Shen 
and  Peterson19  reported  measured  values  of  the  dynamic 
pressure  coefficient  Cp(t).  Fully  wetted,  time  dependent,  ex¬ 
perimental  pressure  distributions  were  compared  with  results 
from  Giesing’s  method27  for  calculating  unsteady  potential 
flow.  Good  correlation  between  the  predictions  and  the  ex¬ 
perimental  measurements  was  obtained  for  both  dynamic 
pressure  amplitude  and  phase  angle  within  the  range  of 
reduced  frequencies  investigated  (K  «  0.23  to  2.30) 

Let  am  and  oIU  denote  the  steady  and  unsteady  incep¬ 
tion  angles,  respectively.  Let  t(K)  be  the  ratio  of  dynamic 
angular  pressure  gradient  (dCp/doL  versus  static  angular 
pressure  gradient  (dCp/do),  at  a  given  location  on  the  foil, 
namely: 

« -(£).  /($?). 

Based  on  Equation  (12),  Shen  and  Peterson  showed  that 
the  unsteady  inception 

angle  for  a  given  reduced  frequency  K  can  be  obtained  from: 


.  .  cos  ♦ 

»IU  -  »0  +  (®l*  -  °o)  — — 


-  oi  sin  4 


\M^r . 


(14) 


The  Shen  and  Peterson19  experimental  results  showed  that 
cavitation  inception  always  initiated  around  x/c  -  0.02.  This 
result  is  in  good  agreement  with  the  assumption  that  cavita¬ 
tion  inception  occurs  at  the  location  of  Cp^.  Due  to  a  phase 
shift  in  the  pressure  distribution,  a  significant  delay  in  incep¬ 
tion  of  leading  edge  cavitation  was  observed  with  the  foil  in 
oscillation.  It  was  concluded19  that  the  influence  of  hydrofoil 
oscillation  on  inception  at  Cp^  of  leading  edge  sheet  cavita¬ 
tion  can  be  reasonably  predicted  by  Equation  (14). 

It  is  remarked  that  full-scale  propellers  are  generally 
associated  with  high  Reynolds  numbers.  The  locations  of 
boundary  layer  transition  and  CPn)ln  may  be  coincided.12  In 
this  case,  the  pressure  fluctuation  term  ACp,  must  be  add¬ 
ed  in  Equation  (12). 


5.2  Cavitation  Inception  at  Transition 

We  now  consider  the  second  case,  when  cavitaton  in¬ 
ception  occurs  in  the  transition  region  of  the  model.  Aside 
from  the  pressure  field,  cavitation  also  requires  a  time  scale 
in  order  for  nuclei  to  grow.  Experiments  on  axisymmetrical 
bodies  by  Arakeri  and  Acosta,9  and  Huang4  showed  that 
cavitation  inception  can  occur  in  the  boundary  layer  transi¬ 
tion  region  when  the  value  of  ( -  Cp,  -  ACp,)  is  smaller  than 
-  Huang  and  Peterson12  in  dealing  with  steady  cavita¬ 
tion  assumed  that 

OiM  *  -  Gj.fj.  -  A(  -pt 


I  ! 


where  om  and  Cpt,  are  model  cavitation  inception  numbers 
and  potential  flow  pressure  coefficients,  respectively.  On  the 
other  hand  cavitation  inception  at  full  scale  was  assumed  to 
occur  when  o*f  »  -Cp,^.  The  difference  in  these  two 
assumptions  is  due  to  time  scale  associated  with  bubble 
growth  and  flow  regime  associated  with  boundary  layer 
pressure  fluctuations. 

Consider  a  previously  discussed  dynamic  run  with  an 
oscillation  amplitude  of  1.57  deg,  frequency  of  5.5  Hz  and 
free-stream  speed  of  9.75  m/s,  corresponding  to  «  2.4  x 
10*  and  k  -  0.432.  This  is  a  typical  value  of  K  for  a  surface 
ship  propeller  if  the  chord  length  at  0.7  radius  is  used  as  the 
characteristic  length.  The  computed  unsteady  potential  flow 
pressure  coefficients  Cp,  +  Cpu(t)  are  shown  in  Figure  23  for 
several  values  of  tut;  see  Table  2.  The  computations  are  based 
on  Giesing’s  non-linear  unsteady  potential  flow  theory.27 
Recall  that  a  »  3.25  +  1.57  sin  tut.  For  purposes  of  com¬ 
parison,  the  dynamic  and  static  pressure  distributions  are 
shown  in  Figure  24.  This  figure  gives  a  -  4.0  degrees  at  tut 
-  30  degrees.  At  the  same  foil-geometric  angle,  due  to  the 
phase  lag,  the  magnitude  of  the  dynamic  pressures  is  notice¬ 
ably  smaller  than  the  static  pressures  near  the  leading  edge. 

This  result  supports  the  experimental  finding1*  that  leading 
edge  sheet  cavitation  is  delayed  with  the  foil  in  oscillation;  see 
Equation  (14).  The  same  trend  has  also  been  reported  for 
aerodynamic  stall.28 

Figure  5c  suggests  that  the  flow  is  in  the  final  stage  of 
transition  at  10  percent  chord,  Pi0  for  R„c  -  2.4  x  10*. 

Figure  15  shows  that  the  pressure  fluctuations  due  to  advanc¬ 
ing  transition  attain  a  maximum  value  at  P10  around  cut  -  30 
degrees.  The  measured  rms  pressure  fluctuations  give  ACpt  - 
0.0145.  Figure  7  shows  that  the  measured  frequentiy- 
occuring  large  negative  pressure  fluctuations  are  approx¬ 
imately  2.7  times  larger  than  the  rms  value  in  the  transition 
region.  If  cavitation  occurs  at  the  natural  transition  point  of 
the  model,  we  have  oyu  -  0.98  *  0.04  -  1.02.  A  much  higher 
value  of  ACpt  -  0.25  may  be  used  if  transition  is  caused  by 
laminar  separation.  This  value  of  o,m  is  smaller  than  the  value 
of  -CPmin.  In  the  previous  cavitation  tests  reported  in 
Reference  19  with  leading  edge  sheet  cavitation,  inception 
always  occurred  near  the  location  of  CPmin.  However,  cavita¬ 
tion  inception  on  axisymmetric  bodies  by  Arakeri  and  Acosta, 
and  Hqapg,  did  occur  in  the  transition  region  where  the  quan¬ 
tity  ( - Cpt,  -  ACpt)  is  still  smaller  than  -  Cp^. 

Further  studies  required  notably  in  three  areas  to 
determine:  (I)  why  the  magnitude  of  pressure  fluctuation 
terms  associated  with  natural  transition  and  laminar  separa¬ 
tion  are  smaller  in  the  present  measurements  than  the  values 
measured  by  Huang  and  Hannan  >«  with  different  in¬ 
struments;  (2)  theoretically  how  the  location  of  transition 
moves  periodically  with  the  foil  in  oscillation;  and  (3)  the 
validity  of  the  assumption  used  above  that  the  location  of 
transition  measured  at  a  steady  mean  foil  angle  is  the  critical 
location  to  trigger  cavitation  inception  when  the  foil  is  in 
oscillation.  This  assumption  requires  further  verification. 

Two  different  types  of  cavitation  inception  phenomena 
have  been  considered  in  this  study.  One  type  is  associated 
with  cavitation  inception  at  the  location  of  CPmin.  The  other 
type  is  associated  with  cavitation  inception  in  the  transition 
region.  Headform  experiments  by  Arakeri  and  Acosta,9  and 
Huang4  were  correlated  with  the  second  type  of  cavitation  in¬ 
ception.  On  the  other  hand,  the  leading-edge  sheet  cavitation 
on  a  hydrofoil  observed  by  the  present  authors  corresponds  to 
the  first  type  of  inception.  It  is  possible  that  a  cambered 
hydrofoil  with  a  smaller  angle-of-attack,  namely  with  a  less 
severe  suction  pressure  peak,  could  encounter  cavitation  in¬ 
ception  at  transition  as  observed  by  Kuiper29  in  his  model 
propeller  experiments. 
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Figure  23  —  Dynamic  Pressure  Distributions  with  a  =  3.25 
+  1.57  sin  cot 
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TABLE  2  -  DYNAMIC  PRESSURE  COEFFICIENT  Cpjt) 


Due  to  the  existence  of  these  two  different  types  of 
cavitation  inception,  Huang  and  Peterson12  computed  a 
significant  scale  effect.  They  provided  a  method  for  correc¬ 
ting  model/full-scale  propeller  cavitation  scaling  in  a  steady 
flow.  The  present  work  is  intended  to  provide  needed  infor¬ 
mation  to  compute  cavitation  scaling  corrections  in  unsteady 
flow.  It  is  shown  that  a  diagram  of  mean  pressure  coefficient, 
Cp,  and  low  frequency  dynamic  pressure  coefficient.  Cp,(t) 


versus  chordwise  location  must  be  constructed  first.  Second 
ly,  the  location  of  boundary  layer  transition  with  the  foil  in 
oscillation  must  be  determined  either  theoretically  or  ex¬ 
perimentally.  Thirdly,  the  magnitude  of  boundary  layer 
pressure  fluctuations  associated  with  frequently  occurring 
negative  pressure  is  determined.  The  selection  of  this 
magnitude  depends  on  the  Reynolds  number  and  amplitude  of 
oscillation  angles  as  shown  in  Figures  8  to  24. 
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Figure  24  —  Comparison  of  Dynamic  and  Statis  Pressure 
Distributions  at  a  »  4.0  Degrees 


6.  CAVITATION  INSTABILITY  AND  NOISE 
6.1  Foil  Oscillation  and  Cloud  Cavitation 

Unsteady  sheet  cavitation  has  been  recently  reviewed 
by  the  15th  ITTC  Cavitation  Committee.  The  emphasis  of  this 
portion  of  the  paper  is  on  the  salient  features  of  sheet  cavity 
instability.  Tanibayashi30  provided  an  insight  into  the  subject 


with  his  description  of  cavitation  on  a  propeller  in  both  uni¬ 
form  and  nonuniform  flow.  He  concluded  that  the  presence  of 
sheet  and  bubble  cavitation  in  nonuniform  flow  can  be 
predicted  by  quasi-steady  methods,  but  that  the  collapse  proc¬ 
ess  cannnot  be  predicted.  Unfortunately,  the  details  of  the 
collapse  process  are  the  controlling  factor  in  erosion,  noise, 
and  induced  structural  vibration.  When  the  sheet  cavity  pro¬ 
duces  "cloud”  cavitation,  erosion,  noise  and  vibration  are 
observed  to  significantly  increase  in  magnitude.  For  example, 
Chiba  and  Hoshimo31  found  that  the  induced  hull  pressure 
had  superimposed  upon  it  pressure  impulses  produced  by  the 
cloud  cavitation  formed  from  the  breakup  of  the  propeller 
sheet  cavity.  A  physical  description  of  cloud  cavitation  and  its 
formation  has  been  given  as  follows:18 

1 .  A  large  portion  of  the  sheet  cavity  becomes  highly 
distorted  and  undergoes  a  significant  increase  of 
overall  cavity  height  in  the  distorted  region. 

2.  Once  this  distorted  region  begins  to  separate  from 
the  main  part  of  the  sheet  cavity,  the  upstream  por¬ 
tion  of  the  sheet  cavity  develops  a  smooth  surface 
and  reduced  thickness. 

3.  The  separated  portion  of  the  sheet  develops  the  ap¬ 
pearance  of  a  cloud  and  moves  downstream  and  ex¬ 
pands  away  from  the  foil  surface.  The  trailing  edge 
of  the  smooth-surfaced  region  then  moves  down¬ 
stream,  becomes  unstable  at  its  trailing  edge,  and 
quickly  develops  the  characteristic  appearance  of 
the  leading  edge  sheet  cavity  elsewhere  along  the 
span.  Alternately,  the  trailing  edge  of  the  smooth 
portion  of  the  sheet  cavity  moves  upstream  to  the 
foil  leading  edge  as  the  cavity  collapses  and  disap¬ 
pears.  Photographs  depicting  this  process  can  be 
found  in  reference  19. 

Ito,32  reporting  one  of  the  first  detailed  experiments 
on  the  subject,  compared  unsteady  cavitation  on  propellers  in 
a  wake  field  with  pitching  three  dimensional  hydrofoils.  He 
concluded  that  the  reduced  frequency  for  the  blade  element  of 
a  propeller  in  a  wake  field  has  an  important  influence  on  the 
formation  of  cloud  cavitation.  The  implication  of  Ito's  work  is 
that  the  wake  field  and  the  propeller  should  be  considered 
together  to  minimize  the  adverse  effects  due  to  the  formation 
of  cloud  cavitation. 

Later  work  by  Miyata,  et  al.16  with  oscillating  two- 
dimensional  hydrofoils  instrumented  with  surface  mounted 
pressure  gages  showed  that  unsteady  wing  theory  was  useful 
in  explaining  the  relationship  between  the  time-dependent 
pressure  distribution  and  cavitation.  They  concluded  that  the 
cavity  collapse  process  is  strongly  influenced  by  the  unsteadi¬ 
ness  of  the  pressure  field  and  the  reduced  frequency 
associated  with  foil  oscillation. 

The  present  authors18  provided  further  details  on  the 
instability  of  sheet  cavitation  and  the  formation  of  cloud 
cavitation.  Their  experiments  were  done  with  an  oscillating 
two  dimensional  hydrofoil  over  a  Reynolds  number  range 
from  1.2  x  10*  to  3.6  x  10*  and  reduced  frequencies  up  to  2.3. 
The  results  indicated  that  the  principal  controlling  param¬ 
eters  were  reduced  frequency,  K,  cavitation  number,  o,  and 
foil  oscillation  amplitude,  oj.  The  maximum  cavity  length, 
lm/c,  is  a  function  of  these  three  parameters  and  cannot  be 
predicted  on  the  basis  of  K  *  0  conditions.  The  role  of  reduc¬ 
ed  frequency  can  be  demonstrated  in  the  following  example. 


For  constant  o,  it  is  possible  to  have  marginal  or  no  cloud 
cavitation  at  some  finite  K,  even  though  it  was  present  at  the 
maximum  unsteady  angle  when  K  -  0. 

The  importance  of  reduced  frequency  has  also  been 
shown  by  Matveyev  and  Gorshkoff.33  They  reported  that  for 
propellers  in  a  uniform  flow  field,  sheet  cavitation  was  less 
noisy  than  bubble  cavitation.  However,  when  the  propellers 
were  in  a  nonuniform  flow  field,  the  sheet  cavities  become 
more  unsteady  and  sheet  cavitation  then  was  noisier  than 
bubble  cavitation.  As  the  work  of  Matveyev  and  Gorshkoff 
points  out,  it  is  of  crucial  importance  that  a  similarity  in 
cavitation  time  history  exist  between  model  and  full  scale. 

The  noise  scaling  relations  assume  this  similarity  as  a  founda¬ 
tion  and  thus  the  importance  of  work  such  as  reported  by 
Bark  and  van  Berlekom1*1  is  reenforced. 

Bark  and  van  Berlekom  tried  to  assess  the  cavitation 
noise  produced  by  a  propeller  in  nonuniform  flow  by  studying 
the  cavitation  noise  generated  by  a  pitching  hydrofoil.  Based 
on  a  correlation  of  photos  of  cavity  life  cycles  and  their 
associated  radiated  noise,  they  concluded  that  good  simulation 
of  specific  events  is  important  and  that  these  important 
events  are  not  generally  described  by  simple  parameters  (e.g. 
collapse  time,  Tc;  variation  of  cavity  area  with  time,  Aft);  and 
maximum  cavity  volume,  Vmax).  One  of  the  most  important 
aspects  of  the  process  is  the  separation  of  the  cavity  (i.e. 
cloud  cavitation  formation)  which  must  be  correctly  scaled. 
They  found  that  cloud  cavitation  formation  can  occur  at  an 
early  stage  of  the  sheet  cavity  collapse.  Furthermore,  in 
agreement  with  the  results  of  the  present  authors19  a  com¬ 
bination  of  long  (chordwise)  cavity  length  and  high  reduced 
frequency  causes  extensive  cloud  cavitation.  As  can  be  con¬ 
cluded  from  the  data  presented  in  their  paper,  large  noise 
transients  were  associated  with  the  cloud  cavitation. 

6.2  Cloud  Cavitation  and  Noise 

The  observed  phenomena  to  be  discussed  in  this  paper 
can  best  be  described  by  reference  to  examples  shown  in 
Figure  25.  Table  3  summarizes  all  of  the  parameters 
associated  with  the  tests  reported  here.  Tolimit  the  scope  of 
the  test  program,  air  content  of  the  water  was  not  varied. 

The  air  content  was  measured  with  70  percent  saturation  in 
reference  to  atmospheric  pressure  at  a  water  temperature  of 
22.2°C  and  tunnel  pressure  of  103.6  KPa.  For  a  velocity  of 
16.4  m/s,  a  =_3.25  +  0.95  sin  cut,  o  =  1.21,  a  plot  of  relative 
sound  power  P2  versus  K  shows  the  existence  of  a  “noise 
bucket."  When  a  =  4.3°  and  K  =  0,  extensive  cloud  cavita¬ 
tion  is  developing  from  the  sheet  cavity  and  the  noise  level  is 
high.  As  the  foil  is  oscillated,  a  leading  edge  sheet  cavity  ex¬ 
periences  an  inception,  growth,  and  collapse  cycle  related  to 
the  impressed  pressure  distribution  of  the  foil.  At  low  reduc¬ 
ed  frequencies  no  cloud  cavitation  is  produced,  the  sheet  cavi¬ 
ty  collapses  toward  the  foil  leading  edge,  and  there  is  a 
significant  reduction  in  the  noise  relative  to  the  condition  of  a 
«  4.3  and  K  «  0.  As  the  reduced  frequency  is  increased  fur¬ 
ther,  a  cloud  cavity  is  produced  during  the  collapse  of  the 
sheet  cavity  and  the  noise  level  again  increases  due  to  cloud 
cavitation  collapses.  This  variation  in  P-  as  a  function  of  foil 
angle  is  shown  in  oscillograph  records.  Figure  26,  for  2  foil 
oscillation  periods  and  in  Figure  27  as  a  mean  noise  variation 
based  on  the  average  of  the  cycles  occurring  in  a  40  second 
period. 

Figurt  25  also  presents  data  that  demonstrate  the 
strong  depend  -nee  of  P2  on  the  water  velocity.  For  o  =  4.3, 

K  »  0,  and  o  =  1.13,  heavy  cloud  cavitation  is  present  and 


the  noise  level  increases  by  a  nominal  factor  of  30  when  the 
velocity  is  increased  from  11.5  m/s  to  16.4  m/s.  Once  foil 
oscillation  starts  the  amount  of  cloud  cavitation  is  significant¬ 
ly  reduced  and  the  velocity  difference  appears  to  have  far  less 
impact  on  the  radiated  noise.  At  large  K  when  the  cloud 
cavitation  is  produced  upon  sheet  cavity  dependence,  the 


Figure  25  —  Influence  of  V*  and  Reduced  Frequency  on 
Relative  Sound  Power  for  a  =  3.25  +  0.95  sin  cot 


strong  influence  of  velocity  is  again  apparent. 

The  influence  of  velocity  on  the  radiated  sound  power 
is  just  as  dramatic  when  the  amplitude  of  oscillation  is  in¬ 
creased  from  0.95  to  1.55  degrees,  keeping  o  constant,  as 
shown  in  Figure  28.  There  it  is  seen  that  with  the  larger 
amplitude  of  oscillation  cloud  cavitation  is  present  over  '  he 
full  range  of  reduced  frequencies.  Although  the  violence 
associated  with  the  cloud  cavitation  at  large  reduced  frequen¬ 
cies  limited  the  ability  to  collect  data,  it  is  apparent  that  P2 
has  a  strong  dependence  on  K  when  the  cloud  cavitation  is 
present.  Reference  to  Figure  25  shows  that  this  latter  effect 
is  not  present  when  little  or  no  cloud  cavitation  is  present. 


TABLE  3  -  TEST  RUN  NUMBERS  AND 
ASSOCIATED  PARAMETERS 

The  variation  in  P2  during  the  foil  oscillation  period  is  shown 
by  oscillograph  records  in  Figure  29  and  as  a  mean  variation 
based  on  the  average  of  the  cycles  occurring  in  a  40  second 
period  in  Figure  30.  Based  on  the  limited  data  presented  in 
Figures  25  to  30,  it  appears  that  for  cloud  cavitation  _ 
originating  from  an  unstable  leading-edge  sheet  cavity  P'2  ~ 
V**  when  o,  K,  on,  oi  are  constant 


The  maximum  cavity  volume,  area  coverage,  etc.,  of  a 
leading  edge  sheet  cavity  have  in  the  past  been  used  as 
parameters  associated  with  the  magnitudes  of  the  cavitation 
noise.  Figure  31  presents  the  maximum  length  achieved  by 
the  leading  edge  sheet  cavity  as  a  function  of  the  reduced  fre¬ 
quency.  From  these  data  and  the  noise  data  of  Figure  25,  it  is 
clear  that  the  maximum  cavity  length  has  essentially  no  cor¬ 
relation  with  P2.  This  basic  conclusion  was  also  deduced  by 
Bark  and  Berlekom.18  As  has  already  been  shown,  the  prin¬ 
cipal  source  of  noise  is  the  collapse  of  the  cloud  cavitation 
generated  by  the  sheet  cavity.  As  shown  in  Figure  32,  for 
small  reduced  frequencies,  the  sheet  cavity  disappears  after 
the  radiated  noise  level  has  peaked.  However,  as  the  reduced 
frequency  is  increased  and  the  formation  of  cloud  cavitation 
is  progressively  delayed  until  it  occurs  at  the  sheet  cavity 
desinent  condition,  the  peak  amplitude  of  the  noise  occurs 
after  the  sheet  cavity  has  disappeared.  In  fact,  for  the  reduc¬ 
ed  frequency  of  1.65,  maximum  radiated  noise  does  not  occur 
until  the  foil  reaches  its  minimum  angle  of  attack.  Before  the 
cloud  cavitation  completely  disappears,  the  inception  of  a  new 
leading  edge  sheet  cavity  has  occurred. 
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Figure  26  —  Cavitation  Noise  Signals  V P2  and  Foil  Angles 
for  Runs  141 1  and  1416,  o  =  3.25  +  0.95  sin  o>t 
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Figure  27  —  Relative  Sound  Power  P2(o)  as  a  Function  of 
Foil  Angle  for  Runs  1411  to  1416  with 
<*,  Im  &  od  for  Sheet  Cavitation 


Figure  26  -  Influence  of  x  and  Reduced  Frequency  on 
Relative  Sound  Power  for  a  *  3.25  +  1.55  sin  uit 
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Figure  30  —  Relative  Sound  Power  P'^tr)  as  a  Function  of  Foil 
Angle  for  Runs  1401  to  1406 


Figure  31  —  Measured  Maximum  Cavity  Length 
vs.  Reduced  Frequency 
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Figure  32  —  Measured  Foil  Angles  at  Maximum  Lift,  Maximum 
Cavity  Length,  Noise  Peak  and  Leading  Edge  Sheet 
Cavity  Desinence  for  Runs  1401  to  1406 
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Figure  33  —  Relative  Sound  Pressure  P2  Over  One  Cycle  for  a 
Continuously  Cavitating  Hydrofoil 
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In  many  instances,  partially  cavitating  hydrofoils  are 
subjected  to  conditions  that  can  be  effectively  simulated  by 
small  amplitude  oscillation  during  which  the  leading  edge 
sheet  cavitation  is  continuously  present.  Based  on  results  for 
the  intermittent  sheet  cavity,  it  is  known  that  sheet  cavity 
stability,  and  hence  propensity  to  produce  cloud  cavitation,  is 
dependent  on  the  reduced  frequency.  The  noise  level  variation 
associated  with  this  type  of  cavitation  is  found  in  Figure  33. 
As  with  intermittent  sheet  cavitation,  the  noise  level  is  low 
when  cloud  cavitation  is  not  present.  Once  cloud  cavitation 
forms,  then  there  is  a  dramatic  increase  in  the  mean  sound 
power  level  and  in  the  time  variation  of  sound  power  level,  as 
shown  in  Figures  33  and  34  respectively.  It  is  apparent  from 
these  results  that  the  desinent  condition  for  the  sheet  cavity 
is  not  required  for  high  radiated  noise  levels,  but  rather  the 
existence  of  conditions  that  promote  the  formation  of  cloud 
cavitation. 

If  a  leading  edge  sheet  cavity  is  considered  as  similar 
to  separated  flow  at  a  foil’s  leading  edge,  then  some  parallels 
can  lie  drawn  with  the  vast  body  of  recently  published  data 
on  dynamic  stall  For  example,  the  following  conclusion  from 
McAlister  and  Carr14  closely  parallels  the  description  of  cavi¬ 
ty  stability  given  by  Shen  and  Peterson. 111 

“The  free-shear  layer  that  was  created  between  the 
region  of  reversed  flow  and  the  inviscid  stream  was  not 
stable  This  instability  resulted  in  a  transformation  of  the 
free  shear  layer  into  a  multitude  of  discrete  clockwise  vor 
tices,  'out  of  which  emergeil  a  dominant  “shear-layer 
vortex  ”  McAlister  and  Carr  go  on  to  further  descritie  the 
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Figure  34  —  Cavitation  Noise  Signals  V P2  and  Foil  Angles 
for  Runs  1510  and  1515,  a  -  3.25  +  0.32  sin  wt 
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upstream  movement  of  a  thin  layer  of  reversed  flow  along  the 
foil  surface.  When  this  reversed  flow  reached  the  leading 
edge,  "...  a  protuberance  appeared  over  the  first  6  percent 
of  the  surface  in  response  to  the  sudden  influx  of  fluid.  This 
protuberance  grew  and  eventually  developed  into  the  'dynamic- 
stall  vortex’  that  has  been  observed  in  high  Reynolds  number 
experiments."  This  dynamic  stall  vortex  moves  downstream 
and  away  from  the  airfoil  surface  just  as  cloud  cavitation  does 
on  a  cavitating  hydrofoil.  Furthermore,  if  the  location  of  flow 
reversal  is  plotted  against  airfoil  angle,34  it  is  seen  that  foil 
oscillation  (i.e.  K  =  0.25)  will  suppress  the  forward  movement 
of  the  reversed  flow  region.  This  is  again  similar  to  the  cavity 
stabilization  at  low  reduced  frequencies  relative  to  the  sta¬ 
tionary  hydrofoil. 

These  types  of  analogies  must  be  used  with  great  care. 
For  example,  one  of  the  conclusions  of  Telionis  and  Koromilas35 
from  their  study  of  unsteady  laminar  separation  is  that 
separation  is  not  affected  by  the  amplitude  of  oscillation.  The 
parallel  with  cavitating  flows  may  break  down  here  due  to, 
among  other  reasons,  the  inertial  considerations  of  growing 
cavities.  Telionis  and  Koromilas  also  have  concluded  that  for 
finite  oscillation  frequencies,  the  point  of  reverse  flow  is 
shifted  downstream  from  the  quasi-steady  location. 

7.  CONCLUDING  REMARKS 


natural  transition  region  than  in  a  fully  established 
turbulent  region. 

With  respect  to  cavitation  noise  the  following  conclu¬ 
sions  were  derived  from  water  tunnel  hydrophone  measure¬ 
ments: 

1.  The  stability  of  a  leading  edge  sheet  cavity  deter¬ 
mines  the  extent  to  which  cloud  cavitation  is  pro¬ 
duced. 

2.  The  formation  of  cloud  cavitation  during  the  life  of 
a  sheet  cavity  is  suppressed  at  small  reduced  fre¬ 
quencies  and  correspondingly  the  noise  level  is  low. 

3.  At  high  reduced  frequencies,  extensive  cloud  cavita¬ 
tion  is  formed  during  the  final  phase  of  sheet  cavity 
collapse  and  the  noise  level  is  significantly  increas¬ 
ed.  Thus  P2  plotted  against  K  shows  a  “noise 
bucket."  When  the  sheet  cavity  is  continuous  with 
time,  a  similar  noise  bucket  is  apparent. 

4.  When  o,  K,  a„,  and  o\  are  kept  constant,  the  in¬ 
fluence  of  velocity  on  js  found  to  be  very  large  if 
cloud  cavitation  is  present.  From  the  limited  data 
obtained  from  these  experiments,  the  sound  power 
associated  with  cloud  cavitation  appears  to  vary  as: 


Depending  on  the  cavitation  resistance  of  the  liquid, 
cavitation  inception  on  a  model  may  occur  either  (1)  at  the 
location  of  Cpmln,  or  (2)  in  the  transition  region  of  the  model. 
In  the  present  experiments  for  a  hydrofoil  with  a  large  suc¬ 
tion  peak,  leading  edge  sheet  cavitation  was  observed  to  take 
place  at  CPmm.  In  this  case,  the  boundary  layer  was  laminar  at 
the  location  of  Cp^  and  an  unsteady  potential  flow  theory 
was  shown  to  provide  a  good  correlation  between  prediction 
and  experimental  measurements  of  cavitation  inception. 

For  the  case  of  cavitation  occurring  in  a  transition 
region,  it  is  shown  that  a  diagram  of  mean  pressure  coeffici¬ 
ent  and  low  frequency  dynamic  pressure  coefficient  versus 
chordwise  location  must  be  constructed  first.  Secondly,  the 
location  of  boundary  layer  transition  with  the  foil  in  oscilla¬ 
tion  must  be  determined  either  theoretically  or  experimental¬ 
ly.  Thirdly,  the  magnitude  of  boundary  layer  pressure  fluctua¬ 
tions  associated  with  frequently  occurring  negative  pressure 
is  determined.  A  general  theory  on  boundary  layer  properties 
with  a  hydrofoil  in  oscillation  is  not  yet  available.  The  present 
work  is  intended  to  provide  some  needed  information  on  this 
subject:  the  location  of  transition  and  the  magnitudes  of 
boundary  layer  pressure  fluctuations.  Experimental  results 
show: 

1.  The  movement  of  boundary  layer  transition  can  be 
detected  by  the  measurement  of  pressure  fluctua¬ 
tions  on  the  foil  surface. 

2.  The  development  of  the  foil's  boundary  layer  is 
delayed  with  an  increase  in  oscillation  frequency. 

3.  The  pressure  peaks  associated  with  advancing  tran¬ 
sition  are  independent  of  oscillation  amplitude  and 
frequency  and  are  identical  with  the  values  meas¬ 
ured  in  the  transition  regions  on  a  stationary  foil. 

4.  In  a  fully  established  turbulent  region,  the  mag¬ 
nitudes  of  pressure  fluctuations  are  independent  of 
oscillation  amplitude  and  frequency  and  identical 
with  the  values  measured  in  the  turbulent  region  on 
stationary  foil. 

5  The  magnitude  of  pressure  fluctuations  measured 
on  the  present  foil  gives  higher  intensity  in  a 


F5  ~  V&o 

5.  The  amplitude  of  oscillation,  o,,  has  an  influence  on 
the  noise  in  that  larger  amplitudes  promote  cloud 
cavitation  formation.  However,  the  limited  data 
available  do  not  permit  more  detailed  discussion  of 
this  point. 

6.  The  cavitation  noise  generated  by  a  stationary  foil 
is  not  indicative  of  the  noise  produced  when  the  foil 
is  oscillated. 

7.  The  reduced  frequency  parameter,  K,  does  not  in¬ 
clude  consideration  of  gross  cavity  dynamics.  Thus 
it  is  not  sufficient  to  predict  the  influence  of  free- 
stream  velocity  on  the  gross  stability  of  the  leading 
edge  sheet  cavity.  This  influence  of  velocity  should 
be  investigated  further  in  order  to  establish  the 
critical  K  at  which  the  radiated  noise  increases  at 
full  scale  speeds. 

The  subject  of  unsteady  cavitation  is  a  topic  just  in  its 
embryonic  state  of  understanding.  Recent  research  confirms 
what  has  been  empirically  observed  for  many  years,  that  is, 
reasonable  simulation  of  erosion,  noise  and  induced  hull  vibra¬ 
tion  requires  a  simulation  of  the  flow  field  in  which  the 
hydrofoil  or  propeller  operates.  A  propeller  operated  in 
uniform  flow  cannot  simulate  some  of  the  critical  details  of 
the  inception,  growth  and  collapse  process  of  leading  edge 
sheet  cavities. 
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